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Using molecular dynamics simulations, we study a liquid model which consists of particles inter- 
acting via a spherically-symmetric two-scale Jagla ramp potential with both repulsive and attractive 
ramps. The Jagla potential displays anomalies similar to those found in liquid water, namely ex- 
pansion upon cooling and an increase of diffusivity upon compression, as well as a liquid-liquid (LL) 
phase transition in the region of the phase diagram accessible to simulations. The LL coexistence 
line, unlike in tetrahedrally-coordinated liquids, has a positive slope, because of the Clapeyron rela- 
tion, corresponding to the fact that the high density phase (HDL) is more ordered than low density 
phase (LDL). When we cool the system at constant pressure above the critical pressure, the hy- 
drodynamic properties rapidly change from those of LDL-like to those of HDL-like upon crossing 
the Widom line. The temperature dependence of the diffusivity also changes rapidly in the vicinity 
of the Widom line, namely the slope of the Arrhenius plot sharply increases upon entering the 
HDL domain. The properties of the glass transition are different in the two phases, suggesting that 
the less ordered phase (LDL) is fragile, while the more ordered phase (HDL) is strong, which is 
consistent with the behavior of tetrahedrally-coordinated liquids such as water silica, silicon, and 
BeF 2 . 



I. INTRODUCTION 

An open question of general interest concerning liquid 
water is the relation between a liquid- liquid (LL) phase 
transition and the dynamic properties Q, 0, IS EL IE 
0, H • The LL phase transition may have a strong effect 
on the dynamic properties of supercooled water, includ- 
ing the glass transition H [^}. In deeply supercooled 
states, some glass-formers show "strong" behavior with 
a well-defined activation e nerg y, while other glass-formers 
display "fragile" behavior [ll| . Water appears to show a 
crossover between fragile behavior at high T to strong be- 
havior at low T0 El El El El The recent stud y on 
the Stillinger- Weber model of silicon E3 , which confirms 
the LL phase transition, suggests that the less ordered 
high-density liquid (HDL) is fragile, while the more or- 
dered low-density liquid (LDL) is strong. These authors 
observed a power-law singularity of the diffusivity in the 
less ordered HDL phase as it approaches the spinodal 
of the LL transition at constant pressure. Recently the 
fragility transition in nano-confined water was studied in 
neutron scatteringexperiments pioneered by the Chen 
group at MIT 0, 3, who found that water appears to 
show a crossover between non- Arrhenius ("fragile") be- 
havior at hig h T to Arrhenius ( "strong" ) behavior at low 
T HH ll2L lilill . Their findings were confirmed using 
NMR by Mallamacc et al. @|. 

A set of realistic water models— ST2, TIP5P, TIP4P, 
TIP3P, and SPC/E — with progressively decreasing tetra- 
hedrality which bracket the behavior of real water have 
been studied to explore the generic mechanisms of 
LL phase transition and anomalies associated with it 
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ity is not a necessary condition for anomalous behavior 
and several spherically-symmetric potentials are indeed 
able to generate density and/or diffusion anomalies. 

In this article, we study a simplified model [Fig. ^ 
based on spherically-symmetric soft-core potentials with 
both attractive and repulsive parts |3l]]. Using exten- 
sive molecular dynamics (MD) simulations, we study 
the static properties, and find that similar to water 
models, the simple Jagla potential exhibits a density 
anomaly (TMD), diffusivity anomaly as well as structural 
anomaly |34| . We further explore the dynamic behavior 
along different constant pressure paths in HDL (more 
ordered) and LDL (less ordered) phase. We observe a 
non- Arrhenius behavior in the less ordered LDL and Ar- 
rhenius behavior in the more ordered HDL phase, as well 
as a dynamic crossover which occurs above the LL crit- 
ical point as the system is cooled down along constant 
pressure paths. We put forth a possible interpretation 
of the dynamic crossover which occurs above the critical 
point in terms of the LL critical point and its associated 
Widom line |5(. Thermodynamic response functions - 
constant pressure specific heat Cp, constant temperature 
compressibility Kt and the structural order parameters, 
Qe (orientational) and t (transnational), support a rela- 
tion between the dynamic crossover and the LL phase 
transition. 

The outline of the paper is as follows. In Sec. II, we 
introduce the spherically-symmetric Jagla potential. In 
Sec. Ill, we describe the method of the MD simulation. 
In Sec. IV, we define the quantities which we study. In 
Sec. V, we investigate the static properties of the model, 
while Sec. VI contains the simulation results of the dy- 
namic properties. Section VII discusses the relation with 
water and another tetrahcdral liquid BcF2. 



2 



II. SPHERICALLY-SYMMETRIC TWO-SCALE 
JAGLA RAMP POTENTIAL 

A spherically-symmetric potential with two different 
length scales has been studied [1 HJ IH HJ HJ |H H3 . 

Here, we study the linear ramp potential but with both 
attractive and repulsive parts |3l| . The potential is de- 
fined as 



U(r) 



oo 

U A + (U A -U R )(r 

U A (c-r)/(c-b) 





b)/(b- a) 



r < a 
a < r < b, 
b < r < c, 
r > c 



(1) 

where Ur = 3.5Uq is the repulsive energy, U a = —Uq is 
the attractive part, a is the hardcore diameter, b — 1.72a 
is the well minimum, and c = 3a is the cutoff at large 
distance [Fig. 




III. METHODS 



FIG. 1: (a) The spherically-symmetric "two-scale" Jagla 
ramp potential with attractive and repulsive ramps. Here 
Ur = 3.5Uq, Ua = —Uo, a is the hard core diameter, 
b = 1.72a is the soft core diameter, and c = 3a is the long 
distance cutoff. 



We apply the discrete MD method |33l |35j, approx- 
imating the continuous potential Eq. Q by step func- 
tions, 



U n {r) 



oo r < a 

U R - k x AU x a + fciAr <r <a+(k x + l)Ar, 
U A b <r < b' , 

U A + kAU 2 b' + (k 2 - l)Ar' < r < b' + k 2 Ar 
r > c 

(2) 

where Ar = 0.02a, Ar' = 0.16a, b' = b + l/2Ar', c = 
c - 1/2 Ar' ', AUi = (U R - U Q )/m with m = 36, AU 2 = 
Uo/n 2 with ri2 = 8, < kx < n\ — 1, and 1 < k 2 < n 2 . 

The standard discrete MD algorithm has been im- 
plemented for part icles interacting with step potentials 
M IM IM El |H Hi|. We use a as the unit of length, 
particle mass m as the units of mass, and Uq as the unit 
of energy. The simulation time is therefore measured in 
units of a-^/m/Uo, temperature in units of Uo/ks, pres- 
sure in units of Uo/a 3 , and density p = Na 3 /L 3 , where L 
is the size of the system and N = 1728 is the number of 
particles. NVE, NVT and NPT ensembles are applied 
in our simulation |2U IS3I • 



IV. THE FOUR QUANTITIES STUDIED 

(a) The diffusion coefficient is defined as 



D = lim 



({r(t> + t)-r(t>)}\ 
6t 



(3) 



where (. . .) t > denotes an average over all particles 
and over all t 1 . 



(b) The static structure factor for wave vector q is 
S(q) = F(q,t = 0), where F(q,t) is the interme- 
diate scattering function defined as 



F(q,t) = (p(q,t)p(-q,0)) 

with the wave vector q = 2iTr/L,i = 1,2, . . 
the Fourier transform of the density function 



(4) 
and 

(5) 



(c) The translational order parameter |3J, |4 



t 



\g(r) - l\dr 



(G) 



where r is the radial distance, g(j-) is the pair corre- 
lation function, and r c = L/2 is the cutoff distance. 
A change in the translational order parameter indi- 
cates a change in the structure of the system. For 
uncorrelated systems, the interaction in the system 
is short-ranged with g(r) = 1, leading to t = 0; for 
long-range correlated systems, the modulations in 
g(r) persist over large distances, causing t to grow. 

(d) The orientational order parameter Q characterizing 
the average local order of the system is defined 
as 



47T x - 

27TT ^ 



\Y t . 



1/2 



(7) 



where Y^ jm (0, (j)) is the spherical harmonic function 
with 9, (j> are angles between the central particle 
and its 12 nearest neighbors, and i ~ 6. In general, 
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the value of Qg increases as the local order of a 
system increases, e.g, Q 6 =0.574 for the fee lattice 
and Qq — 0.289 for uncorrelated systems. 
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FIG. 2: The equation of state P(T,po) of the ramp 
potential for 26 values of po = Na A /L i , where L = 
15.0, 15.2, 15.4, 20.0 is the cell edge. The LL critical point 
(closed circle) is located at P = 0.243, and T = 0.37, well 
above the equilibrium melting line (heavy solid curve). The 
gas- liquid critical point, not shown, is located at much higher 
temperature T gt = 1.446 (P g i = 0.0417 and p gX = 0.102). The 
density anomaly is represented by connecting all the points 
with (dp/dT)p = 0. The locus of the density extrema bounds 
the region of the diffusivity anomaly. 
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FIG. 3: Hysteresis upon compression and decompression 
along a constant temperature path with T = 0.34: (a) Volume 
per particle and (b) Enthalpy per particle. Hysteresis upon 
cooling (C) and heating (H) along constant pressure paths: 
(c) Volume per particle and (d) Enthalpy per particle. 



spinodals. The coexistence line, obtained by the Maxwell 
rule by integrating the isotherms, has a positive slope of 
0.96±0.02/c B a- 3 . 



According to the Clapcyron equation 



V. SIMULATION RESULTS: STATICS 



dP 
dT 



AS 

Ay' 



(9) 



A. Equation of state 

The equation of state of the Jagla model [Fig. is 
obtained using two steps: (1) constant volume simulation 
upon slowly cooling (NVT-ensemble), which allows us to 
obtain the equation of an isochore in a single run; (2) 
constant energy (NVE-ensemble) simulation of individual 
state points of a particular interest. The spinodals are 
defined by the crossing of isochores P(T, p) and P(T, p + 
Ap), where 



dP 
dp 



= 0. 



(8) 



The hysteresis of volume [Fig. 01a)] and enthalpy 
[Fig. |3Tb)] upon compression and decompression along 
constant temperature T = 0.34, and the hysteresis of 
volume [Fig. Etc)] and enthalpy [Fig. Old)] upon heat- 
ing and cooling along constant pressure paths are consis- 
tent with the existence of a LL phase transition. The 
LL critical point, with T c = 0.375, P c = 0.243 and 
p c = 0.37, is located at the maximal temperature on the 



the entropy in the HDL phase is lower than the entropy in 
the LDL phase. Hence, the HDL phase is more ordered 
than the LDL phase, which is the opposite of the LL 
transition found in simulations for water 0] and silicon 
0. 

We also studied the spontaneous liquid-liquid phase 
transition using NPT ensemble [Fig. 0]. Simulations of 
the isobaric heating of the HDL liquid below the critical 
point show that the HDL phase loses its stability and 
spontaneously changes into the LDL in the vicinity of 
the HDL spinodal line [Fig. 0]. The LL critical point 
of the Jagla model, in contrast to those of water and 
silicon models, lies well above the equilibrium melting 
line [Fig. [5] at which the solid and liquid phases coexist 
and are in equilibrium. 

Very recently, Gibson and Wilding |4l| study the fam- 
ily of Jagla potentials with decreasing soft-core distance 
and found that there is a parameter range within which 
the critical point lies below the crystallization line, and 
the coexistence line has a negative slope, resembling the 
situation for water. 
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FIG. 4: Liquid - liquid phase transition. Upon heating along 
constant pressure paths below the critical point, the system 
experience a HDL to LDL phase transition near the HDL 
spinodal line. The inset shows the density anomaly upon 
heating along constant pressure paths. 
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B. Density Anomaly 

The temperature of maximum density (TMD) line is 
defined by (dV/dT)p = 0. Due to the general thermo- 
dynamic relation 



0.15 



8V_ 
df 



dP 
df 



dV 



(10) 



the TMD line coincides with the locus of points satisfying 
(dP/dT)v = 0, which defines the pressure minimum on 
each isochore (Fig. [2J. The density anomaly (density 
increase upon heating along constant pressure paths) can 
also be seen in the inset of Fig. [1] 



C. Thermodynamics 

A sketch of the phase diagram based on the equation 
of state [Fig. |2] is shown in Fig. GJa) and the simulated 
phase diagram is shown in Fig. |5pj) . Different thermo- 
dynamic response functions such as Cp and Kt, which 
diverge at the critical point, have maxima at tempera- 
tures T max (P), which can be regarded as temperatures 
on the extension of the coexistence line above the critical 
temperature (the Widom line la) as we cool the system 
at constant pressure P > P. |42l l43t HI. LHJ . We in- 
vestigate Cp and Kt in the Jagla model along paths a 
and path (3 (Fig. |SJ). We find that, along paths a, Cp 
has a maximum at T max (Po) [Fig. [SJa)], while Kt has 
a maximum at a slightly different T max (Po) [Fig- ED 3 )]. 
The response functions — Cp and Kt — increase continu- 
ously along path (3 [Fig. EI a )] before the system reaches 
the stability limit near the LDL spinodal [la |. 
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FIG. 5: A sketch of the P — T phase diagram for the two-scale 
Jagla model. We study three different paths in the vicinity 
of the LL critical point: (i) Pq > P c (path a). Upon cool- 
ing along path a, the liquid changes from a low density state 
(characterized by a non-Arrhenius dynamic behavior) to a 
high density state (characterized by Arrhenius dynamic be- 
havior) as the path crosses the Widom line, (ii) Pq < P c 
(path 0). Upon cooling along path 0, the liquid remains in 
the LDL phase as long as path (5 does not cross the LDL spin- 
odal line. Thus one does not expect any dramatic change in 
the dynamic behavior along the path (3. (iii) Pq < P c (path 
7). Upon cooling along path 7, the liquid remains in HDL 
phase. Its dynamics, different from that of the LDL phase, 
follows Arrhenius behavior, (b) The P — T phase diagram 
for the two-scale Jagla model. The loci of the specific heat 
maximum Cf? ax , the compressibility maximum A"™ ax , and ori- 
cntional order parameter maximum Q max are similar but not 
quite identical. 



D. Structural order 

Analogously, we can expect that the structural prop- 
erties of the system also change from those resembling 
the LDL phase to those resembling the HDL phase when 
the system crosses the Widom line (the extension of the 
coexistence line). The pair correlation function g(r) for 
different T along a constant pressure path is shown in 
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(b) Compressibility 
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FIG. 6: Response functions for the Jagla model as function 
of temperature for different values of Po > P c [Fig. [S] path 
a], (a) Constant pressure specific heat Cp and (b) constant 
temperature compressibility Kt- Both Cp and Kt have max- 
ima, as is known to occur experimentally |42l [isj . For large 
Po the peaks become less pronounced and shift to higher tem- 
perature as the system is farther away from the LL critical 
point. 
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FIG. 7: (a) The pair correlation function g(r) at constant 
pressure Po = 0.250. The magnitude of the first peak in- 
dicates a HDL-like liquid at low T and LDL-like liquid at 
high T. (b) Distribution of q vectors at constant pressure 
Po = 0.250. The shifts of the first and second peaks in the 
distribution of the q vectors further indicates that the liquid 
changes smoothly from HDL-like to LDL-like as it crosses the 
Widom line, (c) The translational order parameter t and (d) 
the orientational order parameter Qq along constant pressure 
paths. 



rameter Q§ [Fig. Etd)] above the LL critical point along 
path a, indicate that as the system crosses the Widom 
line region [Fig. Efb)], the local structure of the system 
also changes from LDL-like side to HDL-like side, which 
is consistent with the features observed in the thermody- 
namic response functions. 



VI. SIMULATION RESULTS: DYNAMICS 
A. Diffusivity Anomaly 

Figure [8] shows the diffusivity D(p) along seven 
isotherms. There exists a diffusivity anomaly region 
along each isotherm where the diffusivity increases upon 
compressing instead of decreasing. The loci of the dif- 
fusivity extrema where (8D/dP)T = [heavy dashed 
lines in Fig. [21 and Fig. [H] defines the diffusivity anomaly 
region. 
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FIG. 8: Diffusivity as a function of density for seven val- 
ues of temperature. The diffusivity anomaly region, where 
D increases with compression (density), is located within the 
diffusivity extrema lines (heavy dashed line). 

We study the T dependence of D along three different 
constant pressure paths [Fig. EJa)] : 

(i) Path a, Pq > P c (one-phase region); 

(ii) Path /3, Po < Pc (in the LDL phase); and 

(iii) Path 7, P < P c (in the HDL phase). 



Fig. 0a). At low temperature g(r) exhibits a more pro- 
nounced first peak near the hard core distance and a less 
pronounced peak at high T, indicating a change from 
the LDL-like structure to HDL-like structure upon cool- 
ing along paths a. The same behavior can also be seen 
from the static structure factor S(q) [Fig. Efb)]. The 
sharp transition in the translational order parameter t 
[Fig. |7Jc)] and the peaks in the orientational order pa- 



B. Dynamics for P < P c 

For Po < P c , the diffusivity D at high temperature 
T > T c follows the Arrhenius law 

D = D cxp(-^4-) (11) 
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FIG. 9: Dynamic behavior for Jagla potential. The In- 
dependence of diffusivity D along constant pressure paths: 
(a) Po = 0.225 < P c for both path j3 and path 7. The more 
ordered phase (HDL) is strong, while the less ordered phase 
(LDL) is fragile, (b) Path f3 with P = 0.225 < P Cl and 
path a with Po = 0.300 > P c . A dynamic crossover oc- 
curs along constant pressure paths above the critical pressure 
when the Widom line is crossed [Fig. HJ. (c) Path f3 with 
Po = 0.200, P = 0.225 < P c , path 7 with P = 0.225 < P c , 
and path a with P = 0.250, 0.300 > P c [Fig.0. 



with a roughly pressure independent activation energy 
Ea ~ 1.53, while at low temperatures in the two-phase 
region, D behaves differently along path (3 and path 7 
upon cooling at constant pressure [Fig. |SJ. Along path 
(3 [Fig. EJa)] which belongs to the LDL phase, the T 
dependence of D is non- Arrhenius and follows the Vogcl- 



D = D exp — 



B 



Dq exp 



To 1 



T — To J 1 V T-To! , 

(12) 

For P = 0.225 the fitting parameters are B « 0.2, 
To « 0.184, and the fragility parameter / = To/B ps 
0.66 [Fig. E| Quantitatively, the "steepness index of 
the fragility" jig m « 368, with the glass tempera- 
ture T g w 0.192 estimated as the temperature at which 
cxp( J B/(T - T )) = 10 fe , where k = 16 is a typical value 
for thermally-activated systems ^(|. The large value of 
m indicates that the behavior in the LDL phase resembles 
that of a very fragile liquid. Note that the value of the 
index number m largely depends on the Vogel-Fulcher- 
Tamman fitting at the last data point [Fig. EUa)] . 

On the other hand, along paths 7 which belong to the 
HDL phase, D follows Arrhenius behavior with E a w 
6.30, which is much larger than the activation energy at 
high temperature. 



C. Dynamics for P > P c 

For Po > Pc along path a, there is a crossover in 
the behavior of D near T w 0.4 [Figs. E^b) and Etc)]. 
Such a crossover, where Cp passes through a maximum 
[Fig.pfa)!, would be predicted by the Adam-Gibbs equa- 
tion nlEa, 



D = D exp - 



C_ 

TS C 



(13) 



because of the rapid change in the temperature depen- 
dence on S c as the Widom line is crossed. When Cp 
has fallen to a smaller and more slowly changing value, 
the temperature dependence of D assumes an Arrhe- 
nius behavior but with a somewhat larger slope than 
at high temperatures where Fig. HJa) shows Cp to be 
very small. The behavior is similar to what was ob- 
served in experimental studies of the strong liquid BCF2 
and in simulations of S1O2 [l5T |. In fact, the par- 
allel with the case of BeF2 [48( is remarkable. In each 
case the Arrhenius slope extrapolates to an intercept at 
1/T = 0, which is six orders of magnitude above the 
intercept of the high temperature Arrhenius part of the 
plot (which is common to all phases). Thus, the be- 
havior of the HDL-like liquid on the low-temperature 
side of the Widom line can be classified as that of a 
strong liquid. The behavior on the high-temperature side 
of the Widom line, in the LDL-like phase, however, is 
very different, resembling that of the fragile liquid (LDL- 
like), as is clear from Figs. Elb) and Etc). Thus, the 
present spherically-symmetric Jagla ramp potential ex- 
hibits a dynamic crossover from LDL-like (fragile liquid) 
at high-temperature to HDL-like (strong liquid) at low- 
temperature, suggesting the analogous fragile-to-strong 
transition as in water, with the difference that the strong 
liquid is now the HDL phase. 
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VII. DISCUSSION 
A. Jagla Ramp Potential 

The mechanisms underlying the different dynamic be- 
haviors we find can be related to the LL phase transition 
[Fig. 13 a)]. The coexistence line has a positive slope, so 
we have one phase for P > P c and two phases — LDL and 
HDL — for P < P c . According to the Clapeyron equation, 
HDL entropy is lower than LDL entropy, so HDL is more 
ordered than LDL which is the opposite of water. In the 
region of the P-T phase diagram between the LDL and 
HDL spinodals, the system can exist in both the LDL 
and HDL phases, one stable and one metastable. 

(i) The limit of stability of the less-ordered LDL phase 
is determined by the high pressure LDL spinodal 
Pldl(T), which, for our model, is unlikely to be 
crossed by cooling the system at constant pres- 
sure [Figs. Efc) and (d)] — because Plf>l(T) as P c 
for all T except in the immediate vicinity of the 
liquid- liquid critical point [Fig. The dynamic 
behavior of the less ordered LDL phase is non- 
Arrhcnius, which is the characteristic of fragile 
glass formers [Fig. 0. 

(ii) On the other hand, the limit of stability of the more 
ordered HDL phase is determined by the low pres- 
sure HDL spinodal Thdl(-P) [Fig. El , which can be 
crossed by heating the HDL phase at constant pres- 
sure [Fig. ^Jc) and (d)]. That is why the dynamic 
behavior of the more ordered HDL phase can be 
studied only when T < T RDL (P) for P < P c . 

B. Comparison with Water 

For the Jagla model, the dynamic crossover upon cool- 
ing for Pq > P c (Path a) is the same that was observed 
in realistic water models 0, El, which is a fragile to 
strong transition. In water models which have a nega- 
tively sloped coexistence, the density of the high temper- 
ature phase is larger than the density of the low temper- 
ature phase [Fig.^j. The density of the high tempera- 
ture phase is larger than the density of the low temper- 
ature phase if the coexistence line has a negative slope, 
so the high temperature less ordered phase is the HDL 
phase. When water is cooled at Pq > P c [Fig. path 
(3], it crosses the coexistence line but does not cross the 
spinodal of the less ordered HDL phase, which may be- 
come almost horizontal line on the P — T plane as recent 
ST2 simulations suggest Q- Thus, water remains in the 
metastable HDL phase even at very low temperatures, 
since the low-density phase may not nucleate out of the 
HDL phase. Accordingly, above the critical pressure, wa- 
ter may remain fragile even at very low temperatures. In 
contrast, upon cooling at Pq < P c [Fig.ES path a], water 
crosses the Widom line, and its thermodynamic proper- 
ties continuously change from those of the HDL phase 



to those of the LDL phase. This is demonstrated by the 
Cp peak found by cooling water in small pores at atmo- 
spheric pressure [43 ■ Therefore, for P < P c one expects 
a fragile-to-strong transition upon cooling, which is in- 
deed experimentally observed 0, H ■ 




FIG. 10: A sketch of the P — T phase diagram for water 
models, showing patha and /3. Note that unlike the Jagla 
two-scale ramp potential [Fig. ISTa*)] ■ the HDL is less ordered 
than the LDL, so by the Clapeyron relation, the Widom line 
has negative slope. 



C. Comparison with the Tetrahedral Liquid BeF2 

Interestingly, what was observed in a combined 
MD/experimental study of BeF 2 USE^EH show a dy- 
namic crossover similar to what we observed for Jagla 
model for P > P c [Fig. Eta)]. In addition, the MD 
simulations of BeF2 show a density anomaly and spe- 
cific heat maximum close to the point of the dynamic 
crossover at about 2T g . In the Jagla model, the dynamic 
crossover and the Cp maximum occur at higher tem- 
peratures T w 3.5T g . The difference between BcF2 and 
the Jagla model is that a second LL critical point has 
not been directly observed for BeF 2 . Therefore, we can 
not call the region of fast change of the dynamic and 
thermodynamic properties a Widom line. However, the 
extrapolation of the simulation isochorcs in the density 
anomaly region suggests possible existence of a critical 
point at lower temperature and higher pressure. Thus, 
the region of fast changes of the thermodynamic response 
functions is possiblely associated with a Widom line em- 
anating from the hypothetical critical point jSJ E3 hid- 
den in the inaccessible region. As in water, this region 
in BeF 2 has negative slope, suggesting that the dynamic 
crossover in BcF 2 upon cooling is related to the entropy 
decrease from HDL-like on the high-temperature side to 
LDL-like on the low-temperature side. 
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VIII. SUMMARY 

In summary, we systematically study a simple 
spherically-symmetric two-scale Jagla potential with 
both repulsive and attractive parts. We find a LL phase 
transition in an accessible region of the P-T phase dia- 
gram. The Jagla potential also displays water-type ther- 
modynamic and dynamic anomalies, as well as a dynamic 
crossover which occurs as the system crosses the Widom 
line while cooled along constant pressure paths P > Pr . 
Our simulations, similar to simulations of silicon |17j . 
show that the dynamics is Arrhenius in the more ordered 
phase (HDL for Jagla model) and fragile for the less or- 
dered phase (LDL for Jagla model). Our study shows 
that the dynamics is Arrhenius on the low-temperature 
side and fragile on the high-temperature side, as in wa- 
ter. The dynamic crossover for P > P c is consistent with 



(i) the experimental observation in confined geometries 
(small pores) of a fragility transition and (ii) exper- 
imental observation of a peak in the specific heat upon 
cooling water at atmospheric pressure in nanopores [49j • 
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